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A solution based upon small-deflection theory is presented 
for the critical shear stress of curved rectangular panels with 
simply supported edges. Computed curves which cover a wide range 
of panel dimensions are presented; these curves are found to he 
in good agreement with test results. Estimated curves are also 
given for panels with clamped edges . 


INTRODUCTION 


A series of papers has been prepared to provide information 
on the buckling of curved sheet. The problem treated in the 
present paper, which is a part of that series, is the determination 
of the oritical shear stress of a cylindrioally curved rectangular 
panel. For panels having simply supported edges this problem is 
solved theoretically (see appendix), and computed curves are provided 
for finding the oritical shear stress. Estimated results are also 
given for curved rectangular panels having clamped edges. 


RESULTS AND DISCUSSION 


The oritical shear stress t 

cr 

is given by the equation 


for cylindrically curved panels 


cr 


= 


b 2 t 


where 

critical— shear-stress coefficient, established by geometry of 
panel and type of edge support 
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D flexural stiffness of panel per unit length 


E Young's modulus of elasticity 

U Poisson's ratio 

b axial or circumferential dimension of panel, whichever is 

smaller (except where noted) 

t thickness of panel 

Two charts are presented, one for panels having curved sides 
longer than the straight, sides and the other for panels having 
straight sides longer than the curved sides. In eaoh of these 
charts the oritioal— shear-stress coefficient k B is plotted 

against a ourvature parameter Z defined by the equation 

z = “Tv / 1 " H 2 
rt v 

or 



where 

r radius of ourvature of panel 

Panels with simply supported edge s.— The critical-shear-stress 
coefficients k B for curved panels with simply supported edges are 
given in figures 1 and 2. In figure 1, the critloal-shear-etress 
coefficients are given for panels long in the circumferential 
direction; b Is measured axially with the result that the values 
of k a and Z are defined in the manner appropriate to a cylinder. 
In figure 2, the critical— shear-stress coefficients are given for 
panels long in the axial direction; b is measured oiroumferentlally 
so that the values of k and Z are defined in the manner appropri- 
ate to an Infinitely long curved strip. The solid curves In eaoh 
figure are computed, and the dashed curve in figure 2 is estimated 
by a method described in the section entitled "Estimation of Critical 
Stress," 
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With, regard to displacements in the median surfaoe during 
"buckling, the boundary conditions for whioh the curves of figures 1 
and 2 apply are zero displacement along each edge and unrestrained 
motion normal to each edge. The available evidence indicates that 
if both normal and tangential motion were completely restrained, 
the curves of figure 1 would be raised only slightly, but the curves 
of figure 2 might be raised considerably at intermediate values of Z 
for large values of panel length-width. ratio a/b. (See section on 
critical shear stress of ourved panel in reference 1.) 

In figure 1 all the buckling curves for panels of fixed length- 
width ratio approach the curves for a complete cylinder at high 
values of Z. The same trend exists when the axial length is greater 
than the circumferential length, as the curves of figure 2 suggest 
by approaching the slope 0,75 appropriate to cylinder curves at high 
values of Z. This result can be. explained by the fact that 
geometrically a curved panel of any given length-width ratio approaches 
a cylinder as the curvature increases. Panels having a large ratio 
of circumferential length to axial length approach a cylinder at 
lower values of Z (defined as in fig. 1) than panels having a small 
ratio of circumferential length to axial length. 

Estimation of critical stress .— A comparison of the computed 
results for the critical stress coefficients of simply supported 
curved panels (figs. 1 and 2) with results for certain known limiting 
cases indicates the possibility of making reasonable estimates for 
panels of other length-width ratios and edge— support conditions. The 
known limiting cases used in this comparison are the critical stress 
coefficients for the flat panel (Z =0), the complete cylinder, and 
the infinitely long curved strip. (See fig, 3.) Figure h shows the 
comparison for panels long in the circumferential direction. In order 
to permit the comparison, the curve for the strip was plotted by using 
the same parameters as were used in the other curves; that is, k s ard. 
Z were defined in terms of the axial rather than the circumferential 
dimension of the panel (dimension b, in fig. 3). Figure 5 shows 
the comparison for panels long in the axial direction. The cylinder 
ourve is replotted in terms of dimension bg in figure 3 so that 

the same parameters are used as were used for the other curves. 

In each of figures 4 and 5 the first three panel buckling curves 
were computed and the fourth curve was estimated. These estimated” ~ 
curves were obtained by using the known limiting results as guides 
and by extrapolating the trends observed in the cases from which 
computed results were available. 

Panels with clamped edges .— Figures 6 ajid 7 give est ima ted 
theoretical critical— shear-stress coefficients for curved panels 
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with clamped edges. The estimates were made in the aforementioned 
manner by making use of the known shear-stress coefficients for 
cylinders and for long curved strips with damped edges (references 2 
and 3, respectively) and available results for flat panels (refer- 
ences 4 and 5) and by extrapolating from the known results for simply 
supported panels. With respect to edge displacements in the median 
surface during buckling,, the boundary conditions for which the curves 
of figures 5 and 6 apply are zero displacement normal to each edge, 
and unrestrained motion along each edge. The available evidenoe 
indicates that complete restraint of both normal and tangential 
motion would affect the curves of figures 5 and 6 only slightly. 

(See discussion of boundary conditions in reference 1.) 


Experimental Verification 

A study of references 6 to 10, which contain test data on the 
critical shear stress for curved sheet, revealed that in the various 
investigations different types of test specimens (fig. 8) and also 
different methods of defining the experimental oritical stress were 
used. Because of the different types of test specimens and different 
methods used, a wide range of values for the critical stresses must 
be expected. In order to make the comparison between theory and 
experiment, therefore, the tesp data were divided into two groups 
according to whether the buckling load would be appreciably affected 
by initial eccentricities (fig. 9). The theoretical curves used 
were those for simply supported curved panels. 

E xperimental buckling stresses only slightl y affected by in itial 
eccentricities With but three exceptions, the experimental critical 
stresses of Rafel (reference 6) and of Eafel and Sandlin (reference 7) 
correspond to snap buckling (fig. 9(a) )• With any appreciable initial 
eooentx’ioitios, deflections tend to increase gradually with load and 
no snap buckle ocoursj this fact Indicates - that the initial eccen- 
tricities in the test specimens are small. The three specimens In 
which snap buckling did not occur fell in the range of 6<Z<12, 
and the buckling stress was 'taken to correspond to the load at which 
the compressive diagonal stress on one side of the sheet ceased to 
increase with increasing load. The critlaal loads for the Moore and 
Wescoat data were obtained from the loads corresponding to the top 
of the knee of the torque-twist curves ( reference 8) . The torque- 
twist curves represent averages of the behavior of all the pane 3. s in 
the cylinders of reference 8 and are thus relatively insensitive to 
local Imperfections, Eigure 9(a) indicates that buokling stresses 
defined in such a way as to be rather insensitive to local imper- 
fections are in good agreement with, but slightly lower than, the 
theoretical critical stresses for curved panels with simply supported 
edges . 
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E xperi mental Duckling , stresses considerably affected by initial 
eccentricities ".- In the tests "by Chiarito (reference 9) the critical 
stress is defined as the stress at which wrinkling first becomes visu- 
ally perceptible in any panel of a multipanel specimen. Except for 
one snap buckle, the critical stresses in the tests by Kuhn and Levin 
(reference 10) were defined by the point of first departure from a 
straight line of the curve obtained by plotting the readings of a 
diagonally mounted Tuckerman optical strain gage against the applied 
load. In both of these methods initial eccentricities of the test 
specimens can be expected to result in relatively low experimental 
critical stresses. 

A comparison between the theoretical critical stresses for curved 
panels and the experimental critical stresses of Chiarito and of 
Kuhn and Levin is made in figure 9(b) . As might be expected because 
of sensitivity to local imperf actions, the experimental data are, 
in general, considerably below the theoretical curves. The data 

of Kuhn and Levin for = 3 at large values of Z, however, are 

D 

appreciably above the theoretical curves. One possible explanation 
for this rather surprising behavior is that all specimens for ^ » 3 
had two intermediate bulkheads which were not attached to the sheet 
but may have provided additional restraint against buckling at high 
curvatures by preventing inward displacement of the sheet on buckling. 


CONCLUDING- REMARKS 


The critical shear stresses given by a theoretical solution 
based on small-deflection theory for simply supported curved rec- 
tangular panels were found to be in good agreement with experimental 
critical stresses defined in such a way as to be rather insensitive 
to local Imperfections. A method is suggested for estimating the 
critical shear stresses for curved panels having length-width ratios 
and types of edge support different from those for which computed 
results are available. 


Langley Memorial Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., March 11, 19^-7 
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APFEWDIX 

THEOESHICAL SOU3TION 
Symbols 


b 

m, n, p, <i 
r 
t 
u 


w 

X 

y 

D 


E 

F 


axial or circumferential dimension of panel, 
whichever is larger 

axial or circimferential dimension of panel, 
whichever is smaller 

integers 

radius of curvature of panel 

thickness of panel 

displacement of points on median surface of panel 
in axial (x-) direction 

displacement of points on median Burface of panel 
in circumferential (y-) direction 

displacement of points on median surface in radial 
direction; positive outward 

axial coordinate of panel 

circumferential coordinate of panel 


flexural stiffness of panel per unit 
.. >3 

length 


Et~ 


^ 12(1 - ^ 

Young's modulus of elasticity 

Airy's stress function for median- surface stresces 

' d : 


produced by buckle deformation 


shear 


\ 

stress; 2-*-, compressive stress in y-direction; 

& ^ ^ 

— of compressive stress in x-direction 1 

Sy^ / 
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(j|.A - , 2 ) 


s inn 

k a 


curvature parameter 
coefficients of terms in deflection functions 


critical shear-stress coefficient appearing in 
formula t. 




“k- 




e » - 
*b 


3203k 


320 3 k B 


cr 


(p 2 + a 2 0 2 ) 2 + ■ 


h 2 t 


120 ^Z 2 


(p2p2 + g,2)2 + 


it^p 2 + q. 2 0 2 ) 2 
120&p^Z 2 

it lf '(p 2 0 2 + q, 2 ) 2 




Poisson's ratio 



critical shear stress 

s^ + si 

Sz 2 Sy 2 

Si. + 2 ^ - + £. 

Sz^ Sx 2 Sy 2 Sy^ 

* lnverso of defined "by V ^(v^w) 


v 


Method of Solution 

Equation of equilibrium . - The critical shear stress which causes 
a curved rectangular panel to buckle may be obtained by solving the 
following equation of equilibrium (reference l) • 


4 

D V w 


Et -4 c£w + 2 t t - 
r 2 Sz^ cr Sz Sy 


0 


(1) 


where z and y are* respectively } the azial and circumferential 
coordinates. Pi vision of equation (l) by D gives 


( 2 ) 
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The equation of equilibrium may Tbs represented "by . 

Qw ™ 0 

where Q is the oporator dofined hy 


Q “ + 


Ig£ v -4 

b^ bx h 


+ SCg 


£ S 2 

t)2 "cby 


(3) 

(4) 


Solution for panels haying axial., length aroater than circumfer - 
ential length . - Equation (3) may ho solved hy tho Galerkin method as 
outlined in references 1 and 11. As suggested in reference 1 for 
simply supported curved rectangular panels, the following series 
expansion is used for w: 


w «= 


co 00 

y~ r 


a fftn SJSs 


(5) 


The coordinate system used is shown in figure 10(a) . The coefficients 
s w, are then choson to satisfy the equations 



sin £2£ gin 932L Qw dx dy « 0 
a h 


(6) 


where 


p => 1,2, . . . and q « 1,2, . . . 


Whon the operations indicated in equation (6) arc performed, a Bet 
of homogeneous linear algebraic equations in the unknown coefficients 
Bjq is obtained with k B appearing as a parameter. The solution 
for the critical-shear-stress coofficient k s is thon found as the 
^minimum value of k s for which the algohraic equations have a 
nonvanishing solution for the unknowns a^. 

The boundary conditions at each edge that are implied by this 
method of solution aro gero radial deflection and edgo moment, no 
displacement along tho odge, and free displacement normal to tho 
edge (seo referenco 12); that is. 
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■w(o,y) = w{a,y) « 0 


T (0,7) 


^w 


a * 2 “ a * 2 

v(0,y) = v( B.,j) m 0 


(a,y) = 0 


w(x,0) *» v(z,"b) *» 0 

0 fee, - 0 

u(x,0) = u(x,b) = 0 


>(7) 


0 (°,y) - 0 M - 0 


0 - g§ <*»»> 


Substitution of the -values of Q and w given by equations (4) 
and ( 5 ), respectively, into equation (6) leads to the following set 
of algobraic equations: 


■®KL 


(p 2 + q 2 p 2 f + 


V*z 2 




Ap 2 + 9.% 2 ) 2 


®an 


mn-pq 


n=sl (m 2 - p 2 ) (n 2 - q 2 ) 


«= 0 (8) 


-where p = 1,2, . . . ; q «= 1,2, . . . ) and the summation includes 
only those valueB of m and n for which m t p and n t q are 
odd. The condition for a nonvanishing solution of these equations 
is the vanishing of the determinant of the coefficients of the 
unknowns a^. This infinite determinant can be factored into the 

product of two infinite .subdeterminants, one in which p ± q is 
even, and one in which p t q is odd. The vanishing of these 
subdeterminants leads to the following determinants! equations. 

The equation in which p t q is even is 
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?U 

a a 

13 22 

a 

31 

a l5 "24 

a 

33 

a 

17 

a 

26 

•a 

35 • • 

p=l, q.=l 

«u 

° t 

0 

0 4 

0 

0 

4 

35 

0 . . 

P=lj 4=3 

0 

Mno 4 
5 

0 

0 f 

0 

0 

4 

9 

0 . . 

P=2, 4=2 

4 

9 

4 M 
“22 

5 

4 

5 

-2Q 0 
.63 ° 

& 

25 

. 28 
135 

0 

k . . 
7 

P=3; 4=1 

0 

0 

5 

% 

25 

0 

0 

-ji 

175 

0 . . 

P=l, 1=5 

0 

° '% 

0 

M-i c 4o 
'27 

0 

0 

11 

0 . . 

p=2, 4=4 

8 

45 

| 0 

8 

25 

-is ^24 
27 

-22 

35 

.56 

99 

0 

8 

• • 

3 

P=3, 4=3 

0 

0 ^ 

25 

0 

0 

35 

M33 

0 

. 108 
125 

0 . . 

p=l, 4=7 

0 

0--26- 

135 

0 

0 -5$ 
99 

0 

m 17 

*28 

13 

0 . . 

p=2, 4=6 

4 

35 

1 0 

16 

175 

20 o 
11 

108 

125 

28 

13 

M £6 

11 * * 

P=3, 1=5 

0 

0 7 

0 

0 - 
3 

0 

0 

>36 

11 

35 • • 



= 0 (9) 
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and the equation in which p i q. 

is odd is 






®12 

S 21 

^4 

*23 ^2 

a l6 ®25 

a 34 

*18 


P=l> 

1=2 

1 %> 

4 

9 

■ 0 

f 0 

0 1 

0 

0 

28 

135 

P= 2 , 

1=1 

! 1). 
i 9 

“21 

8 

45 

o it 
5 

-55 0 

8 

25 

16 

189 

0 

P=l, 

1=4 

0 

8 

45 

** 

-f ° 

o' g 

0 

0 

56 

99 

P=2, 

<1=3 

1 

0 

8 

7 

^2;? _36 

25 

-1 0 

22 

35 

-2=6 

55 

0 

P=3> 

1=2 

i 

0 

4 

5 

0 


o -f 

0 

0 

28 

75 

P=l, 

1=6 

0 

4 

35 

0 

"I ° 

“lS -§2 
11 

0 

0 

28 

13 

P=2, 

1=5 

go 

63 

0 

4o 

27 

0 -? 

20 Mge 

'll ■ 

8 

3 

80 

117 

0 

P-3, 

1=4 

1 

0 

8 

25 

0 

22 o 

35 ° 

0 *3 

34 

0 

_5§ 

55 

P=l, 

1=8 ! 

0 

16 

189 

0 

-is o 
55 

0 Sf 

0 

“l8 

112 

45 

p=2, 

1=7 

28 

135 

0 

5i 

99 

0 -28 
75 

o 

13 

5i 

55 

-ii 2 

45 

n^rj 


where. 



j (p 2 + i 2 (3 2 ) s 


l2$\ h Z 2 
n\p 2 + 1 2 {3 2 ) 2 
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These determinants give the buckling loads of curved panels with 
various length-width radios (p ^ 1} for buckle patterns respectively 
symmetrical and antisyrmetrical about the center of the panel. 

By use of a finite determinant including the rows and columns 
corresponding to the most important terms in the expansion for w 
(equation (5)), equations (9) and (10) were solved by a matrix 
iteration method (reference 13 ) for the lowest value of k s which 
satisfied those equations* The lower of the two values of k B 
found by solving equations ( 9 ) and (10) is the crltical-shoar-stress 
coefficient for the particular values of (3 and Z under consider- 
ation* The curve of critical-shear-stress coefficient against Z 
for a given value of 3 computed in this manner shews cusps; however, 
as the precise location of all the cusps involves a prohibitive amount 
of labor without any significant increase in accuracy, the cusps were 
faired out in figure 2. Table 1 presents the relative magnitudes of 
the coefficients of the terms used in the solutions, and tablo 2 
gives tho computed stress coefficients* 

Solution for -panels having circumferential length greater than 
axial length . - When the circumferential dimension is greater than the 
axial dimension, a and b can be interchanged in oquation ( 5 ) in 
order to retain b as the shorter dimension, as follows: 

w = ^ r sin sin (ll) 

m=l n-1 ° a 


The coordinate system used is shown in figure 10(b) * This problem 
is solved in a manner similar to that used in solving the problem 
Involving axial dimension greater than circumferential dimension* 
The set of algebraic eq.uations for the unknown Fourier coefficients 
is new 


(p 2 P 2 + q . 2 ) 2 + 


lgpQp^Z 2 


^(p 2 P 2 + <l 2 ) 2 

_B2- ^_S2 . 


+ 3 2 p 3 * s y ~ 

rt 2 m=I n=i 02111 (m 2 - p 2 )(n 2 - q. 2 ) 


Mm. 


( 12 ) 

* 0 


where m t p and n i q. must bo odd for to have a value* The 

condition for a nonvanishing solution of those equations is the 
vanishing of the determinant of the coefficients of tho unknowns a^* 
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This determinant can "be factored into two infinite subdeterminants 
which are identical to those in equations (9) and (10) except for 
the diagonal terms. In this case, each diagonal term is 

replaced by 




(p?-p2 + *2)2 


+ l2 (3 8 p 4 Z 2 

Jt4(p2p2 + q2}2 


The conditions for the vanishing of these determinants were 
found In the same manner as that reed In solving equations (9) 
and (lo) for buckling loads of panels with various length-width 
ratios for buckle patterns symmetrical or antisymmetrical about the 
center of the panel. The lower of the two values of k ? is the 
critical-shear-stress coefficient for the particular values of 3 
and Z under consideration. The relative magnitudes of the coef- 
ficients of the terms used and the computed stress coefficients are 
presented in tables 1 and 2, respectively* Figure 1 shows the 
critical-shear-stress coefficients for simply supported curved 
panels having the circumferential length greater than the axial 
length j the curves were faired in a manner similar to that for the 
curves of figure 2. 
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TABLE 2.- THEORETICAL CRITICAL-SHEAR-STKESS COEFFICIENTS 


k- a — »f . 

ai 


a 

Z 

k 

s 

Even. dete^i 

Odd. determinant 


1 

0 

9-35 



1 

9-44 

11.59 


10 

11.65 

12-77 


30 

18.57 

17-59 


100 

34-65 

33*55 


1000 

157.1* 

164-5 

1.5 

0 

7-07 

7-97 


1 

7.12 

8.03 


10 

8.55 

9.75 


30 

14.30 

15-38 


100 

30. 54 

27-15 


1000 

136.6 

129.7 

2 

1 

6.62 

6.65 


10 

7-65 

8.43 


30 

32.48 

14.29 


loo 

26-96 

26-19 


1000 

117.3 

H8.9 





a 

Z 

k 

s 

Even determinant 

Odd determinant 

1-5 

1 

7.37 

7.99 


10 

IO .38 

9-49 


30 

15-23 

15.51 


100 

32.24 

30.73 

2 

1 

6.68 

6.61 


10 

8.98 

8.95 
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Figure 2.- Critical-shear -stress coefficients for simply supported 
curved panels having axial length greater than circumferential 
length. 
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Figure 3.- Curved rectangular panel. Limiting cases for a complete 
cylinder and an infinitely long strip also shown. 
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Figure 4.- Critical-shear-stress coefficients for simply supported’ 
curved panels compared with those for cylinders and long curved 
strips (k s and Z are defined in terms of axial length of panel) . 
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Fig. 6 



Figure 6.- Estimated theoretical critical-shear-stress coefficients for 
curved panels with clamped edges and having circumferential length 
greater than axial length. 
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Figure 7.- Estimated theoretical critical-shear-stress coefficients for 
curved panels with clamped edges and having axial length greater 
than circumferential length. 
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Fig. 8 
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Figure 8.- Transverse sectional views of specimens used to determine 
critical shear stress of curved panels in various investigations. 
Different length-width ratios of panels were obtained by varying 
bulkhead or ring spacing and in some cases by spacing of 
longitudinal stiffeners. 
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(a) Specimens slightly affected by initial eccentricities. 

Figure 9.- Comparison of theoretical critical-shear-stress coefficients 
of simply supported curved panels with test results of other 
investigations. 
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(a) Panel long in axial direction. 


(b) Panel long in circumferential direction. 


Figure 10.- Coord inat e systems used in theoretical solution. 
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